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Abstract 

Random sequential adsorption is an irreversible surface deposition of ex- 
tended objects. In systems with continuous degrees of freedom coverage follows 
a power law, 9{t) ~ 9j — ct~", where the exponent a depends on the geometric 
shape (symmetry) of the objects. Lattice models give typically exponential 
saturation to jamming coverage. We discuss how such function 9{t) can be 
computed by series expansions and analyzed with Fade approximations. We 
consider the applications of efficient Monte Carlo computer simulation method 
(event-driven method) to random sequential adsorptions with high precision 
and at very long-time scale. 
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1. Introduction 

Random sequential adsorption^ (RSA) refers to a type of irreversible adsorptions 
of spatially extended objects. Two features characterize RSA: (1) the objects, once on 
the surface, are permanently stick to the surface without further thermal relaxation; 
(2) only a single layer is deposited, and the objects have hardcore volume exclusion. 

A standard example^ of RSA is the adsorption of spherical macromolecular parti- 
cles in solution on a glass surface. Some surface reaction dynamics may be described by 
RSA on lattices. In comparison with a corresponding equilibrium system of particles, 
RSA shows unusual properties, most notably the existence of a jamming coverage and 
a super-exponential correlation of particle density. 

One of the central theoretical problems in RSA is the computation of the cov- 
erage 9{t). Such problems are intractable for exact solutions except perhaps some 
one-dimensional systems. ^'^ Series expansion and Monte Carlo simulation offer semi- 
numerical or numerical results which can be improved systematically. They are perhaps 
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the best methods available. In this article, we describe some of the advanced techniques 
of attaining that goal. 

2. Series Expansions 

Earlier calculations^'^ in RSA were based on truncation of an infinite hierarchy of 
rate equations together with mean-field type approximations (and shielding property) 
to close the equations. It is only recently that these equations are systematically used, 
resulting very long series. ^"""^^ 

To illustrate the series expansion method,^ let us consider the deposition of dimer 
on square lattice. The basic quantity of interest in computing a series is the probability 
P{A) that the given set A of sites are empty; we do not care other sites being occupied or 
empty. It is a marginal probability of the full probability distribution. Such probability 
can only decrease. The rate of decrease is proportional to the probability that the 
current configuration A can be destroyed by depositing a dimer with at least one site 
in A. Thus we have 

dPiA) 



oc 



dt 

ways of destroying A 

where the summation is over all possible ways of depositing a dimer at a pair of empty 
sites such that at least one site is in A; A' = A if two of the dimer sites of a deposition 
attempt are within A, or A' is one site more than A so that the deposition with one 
site in A and one site outside A can be carried out. The proportionality constant sets 
the time scale. Without loss of generality, we can take it to be unity. The first two 
equations look like these: 

dP{o) 



dt 



4P(oo), (2a) 



dPioo) ^ p^^^^ + 2P(ooo) +4P(oo). (26) 



dt 

There are four ways to destroy a single empty site, provided that the nearest neighbor 
site is also empty. Using the assumption that initial conditions are lattice symmetry 
invariant, we can write them simply as 4P(oo). The second equation is derived similarly. 

For a general discussion, we write the rate equations symbolically as 

'-^ = jCP{A), (3) 
where £ is a linear operator defined by 

£P{A)^J2'^A'P{A'). (4) 



The n-th derivative is then 
d^'PiA) 



df^ 



£"P(A), with P{A\^^^ = 1 foraUyl'. (5) 
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The initial conditions are such that all sites are empty at t = 0. 



The above scheme can be implemented on computer as follows: starting from some 
initial pattern (or called configuration) we apply £ to get Ap from which we obtain 
the first derivative at t = 0, as dP{A'^)/dt\t=o — cjP(Aj) |t=o = Cj- The second 
derivative is obtained recursively by applying C again to the new patterns as 



The superscript 2 in c| refers to the second generation of patterns. The relationship of 
the patterns and the rate equations can be represented as a tree with arbitrary number 
of branches. Each node of a tree contains a pattern A, and its associated derivatives 
d^P{A)/dt'^ for i = 0, 1, 2, . . . , n. Each node, if its derivative has been computed, has 
a list of pointers to the children nodes and associated coefficients Cj. Thus a node 
together with its children symbolically represents one rate equation. 

If the tree nodes are traversed in a depth-first manner, then we need not store all 
the nodes at the same time. The memory requirement is proportional to the number of 
levels of the tree. Simple counting method usually uses this strategy. The same node 
may have to be expanded many times, and repeated work has to be done. However, 
minimum book-keeping is needed. 

Computation can be performed faster if we do not repeat the same expansion and 
calculation. However, it is necessary then that all nodes are kept; the tree can be 
generated in a breadth-first fashion one level at a time. Equivalent nodes due to lattice 
symmetry are treated as the same node. If a particular node is already generated before, 
a pointer reference is made to the existing node. Such strategy is known as dynamic 
programming in computer science. 

The most practical and efficient algorithm is a combination of several methods. To 
conserve memory, depth-first expansion after every D levels is made; while the highest 
levels use a simple counting algorithm which uses a negligible amount of memory. 

The above method is applicable to a wide range of problems, such as diffusion- 
reaction models, Ising relaxation dynamics, and other dynamic processes described by 
a hierarchy of rate equations. 

Summary of Results 

Series^° for the dimer on square lattice (17 orders), monomer with nearest neigh- 
bor exclusion on square lattice (20 orders), dimer on honeycomb lattice (20 orders), 
and monomer with nearest neighbor exclusion on honeycomb lattice (23 orders) are 
obtained. To our knowledge, there series are the longest known. 

The discrete counting problem for RSA on lattice can be generalized to RSA on 
continuum, where summations are replaced by integrals. We note that these integrals 




(6a) 





k 
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can be related to the cluster integrals with a diagrammatic rule similar to the Mayer 
theory for fluid. A five-order series for disks and a seven-order series for aligned squares 
or cubes are availableJ'^^ 

Pade Analysis of Series 

Given a finite number of terms in a series, how can we say anything about the 
whole time domain t from to oo? Clearly it is not possible in general. Our experience 
with RSA series suggests that we can give reliable estimates from the series for all 
times. ^ 

The method of Pade approximants^^'^'^ is very useful in this respect. Given a series 
f{x) to order L, we determine two polynomials Pn{x) and Qd{x) of degree N and D 
respectively, such that 

^^^)-^ = ^(^'^^')' N + D<L. (7) 
Qd{x) 

It is a powerful way of extending the domain of convergence of the original series. To 
accelerate the convergence further, new variables are introduced,^ e.g., 

s = l-exp(-fe(l-e-*)). (8) 

The functional form is chosen in such a way so that the series in the new variable s 
is most closely resemble the asymptotic behavior of the coverage 9{t) at large t. The 
above specific form is encouraged by the exact solution of the one-dimensional dimer 
problem. In fact we have 6 = s with 6 = 2 in such case. The convergence among various 
Pade approximants is improved greatly by the transformation. 

Here is an example of Pade approximant for the dimer on square lattice 

(^2.962963 s + 0.03206897 - 2.195246 - 1.073721 + 0.9207869 + 0.5556586 

- 0.04386743 - 0.05303456 s^^ / (l + 1-733045 s - 0.2568919 - 1.942572 

- 0.5852424 + 0.7908992 + 0.4421557 - 0.0493306 s'^ - 0.0513337 s^^ , (9) 

where s is given by Eq. (8) with b = 1.35. The result is stable against variation in b. 
Error can be estimated from the convergence of various Pade approximants and Monte 
Carlo simulation results. The above Pade approximant is accurate to 10~^ for alH > 0. 
Such a high degree of accuracy embedded in a simple formula is remarkable. 

3. Monte Carlo Simulation 

Monte Carlo simulation^^"^^ is a simple and useful method to get a quick result 
on an otherwise difficult problem by analytical means. RSA model can be simulated 
rather easily with a simple program. Our emphasis here is how to make Monte Carlo 
simulation more efficient. During late stage of simulation, most of the sites are already 
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occupied or blocked and few adsorption attempts may succeed. The dynamics thus 
becomes very slow, especially for deposition on continuum. 

To overcome such slow dynamics, event-driven algorithms^^"^^'^^ are devised where 
the deposition attempts can be chosen in such a way so that it is accepted with a high 
probability. 

To illustrate the method, let us look at the dimer RSA on square lattice again. At 
late stage, the number of places where a dimer can be placed is small. We can afford 
to enumerate them and store them in a list. In the next stage of simulation, we pick 
a possible dimer position from the list randomly, and make the deposition attempt. 
The difference between the new method and the original is that time step has to be 
different. In particular, each deposition attempt taken from the list should advance the 
time step stochastically by 

= \ , ^"^^ . + 1 ) , (10) 
Ntot VLln(l -r)J y ' ^ ' 

where Ntot is total number of lattice sites, and r is the ratio of the number of ways 
of possible depositions on the list to the total number of ways of depositions. a 
uniformly distributed random number between and 1. The original method has a 
constant time step 1 /A^tot 

Similar idea is used for RSA of disks on continuum,-^* where the area is divided into 
small squares. Each square is classified as available or unavailable for further deposition. 
The classification is nontrivial, but can be done efficiently in two dimensions. The new 
algorithm enable us to simulate the system until jamming state in finite computer time. 
We confirmed the Feder's law^^ 

0{t)-djo,t-^'\ (11) 

and also obtained a very precise value for the jamming coverage 6j ~ 0.547069. 

An important application of the event-driven algorithms in RSA is the deposition 
of more complex objects^^ like polymers and macromolecules represented as random 
walks. A lattice RSA model that we have proposed describes an irreversible adsorption 
of polymer chains in good solvent (self-avoiding random walks). The coverage is a 
very slow function of time. Thus a straightforward simple algorithm can only probe the 
short-time behavior. The event-driven algorithm requires to list all the available ways of 
deposit polymers on the lattice. Clearly such enumeration is not possible for long chains 
due to the huge number of possible configurations a polymer can take. Thus we only 
enumerate the first few steps of the chain, and sample remaining sites stochastically. 
The length of the stored segments of the chains is dictated by the computer memory 
available. The event-driven method substantially extends the time scale, and is able to 
simulate the process over 12 decades in time. Our numerical results suggest that the 
coverage as a function of time obeys 

e-ej(x (12) 
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in a broad time domain for fong chains, where N is the length of the self-avoiding walks. 
The jamming coverage is consistent with 9j oc N"^-^. 

When the objects in RSA deposition are allowed to relax, specifically, a random 
diffusion on the surface, the picture changes a lot. In particular, a full coverage can be 
reached (depends on the geometry of the objects). RSA with diffusional relaxation is 
studied by Monte Carlo simulations^*^ '^^ as well as series expansions. Unlike simple 
RSA, the late stage dynamics typically obeys a power law even on discrete lattices. 
Some of the peculiar behavior^-*^ have not been understood very well. Other research 
direction is the adsorption involving multi-layers. ^^'^^ It is interesting to investigate 
new phenomena and to develop efficient methods for such more complex problems. 
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